Travelling waves in pipe flow 
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A family of three-dimensional travelling waves for flow through a pipe of circular cross section is 
identified. The travelling waves are dominated by pairs of downstream vortices and streaks. They 
originate in saddle-node bifurcations at Reynolds numbers as low as 1250. All states are immediately 
unstable. Their dynamical significance is that they provide a skeleton for the formation of a chaotic 
saddle that can explain the intermittent transition to turbulence and the sensitive dependence on 
initial conditions in this shear flow. 

PACS numbers: 47.20.Ft ,47.20.Lz ,47.35. +i 



Based on decades of studies it is consensus that Hagen- 
Poiseuille flow through a pipe of circular cross section 
belongs to the class of shear flows that does not be- 
come linearly unstable. Nevertheless, it undergoes an 
intermittent transition to turbulence for sufficiently high 
Reynolds numbers and sufficiently large initial pertur- 
bations, as first documented in the classic experiments 
by Reynolds Since then many studies have ana- 

lyzed the mechanisms of transition and the properties 
of the turbulent state H H SH H. Particularly rele- 
vant to the present analysis are the investigations by Dar- 
byshire and MuUin j3] which clearly show a strong sensi- 
tivity to perturbations and a broad intermittent range 
of decaying and turbulent perturbations in an ampli- 
tude vs. Reynolds number plane. Experimental and nu- 
merical studies 0, 1^ show that the turbulent flow in 
the transition region is dominated by downstream vor- 
tices and streaks. Various models for their dynamics 
have been analyzed |0| . Taking the full nonlinearity 
into account WalefFe developed the concept of a nonlin- 
ear turbulence cycle for the regeneration of vortices and 
streaks In addition, stationary solutions and travel- 
ling waves have been found in the full nonlinear equations 
for plane Couette, Tay lor- Couette and plane Poiseuille 
flow 12, 1 3,14", Iff, Tef. It has been suggested that these 
structures provide a skeleton for the transition to tur- 
bulence and the observed intermittency 0, 0| . They 
clearly dominate various observables in low Reynolds 
number turbulent flows jlOj . and are also relevant for 
an understanding of the effects of non-Newtonian addi- 
tives [13. 

The existence of exact coherent states in pipe flow has 
been an object of speculation for some time 5, ^0- As 
we will show here Hagen-Poiseuille flow supports fami- 
lies of travelling waves with structures similar to those 
observed in other shear flows as well. This underlines 
the significance of vortex-streak interactions also in this 
system and opens alternative routes to modelling and 
controlling pipe flow. 
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The existence of stationary states in plane Couette flow 
is connected with an inversion symmetry in the laminar 
profile. In the absence of such a symmetry in pipe flow 
the simplest states we can expect are travelling waves 
(TWs), i.e. coherent structures that move with constant 
wave speed. The wave speed depends on shape and struc- 
ture and is not known in advance. We search for TWs us- 
ing a Newton- Raphson method applied to the full Navier- 
Stokes equations for an incompressible fluid, 

dtu. + (u ■ V)u = I'Au — Vp 
V-u = 0. 

The kinematic viscosity v is kept as a free parameter for 
the continuation method and the Reynolds number is de- 
termined by the mean flow velocity J7, i.e. Re = 2RU /f, 
where R is the pipe radius. For the numerical represen- 
tation of the Navier-Stokes equation we use a Fourier- 
Legendre collocation method in cylindrical coordinates 
with Lagrange multipliers to account for no-slip bound- 
ary conditions at the wall and the constraints that the 
flow field is solenoidal, analytical and regular in the cen- 
ter for r — 0. The program was verified by reproducing 
literature values for the linearized problem (2fli, for the 
nonlinear dynamics of optimal modes |2lj| and for the sta- 
tistical properties of fully developed turbulent flow up to 
Reynolds numbers of 5000 The Newton- Raphson 

method was implemented with a spatial resolution of 
|n/ll| -I- |m/ll| < 1, where n and m are the azimuthal 
and downstream wavenumbers, respectively, and 56 Leg- 
endre polynomials radially. This gives us about 5600 
dynamically active degrees of freedom. 

Initial conditions for the Newton method are obtained 
by embedding pipe flow into a larger, parameter depen- 
dent family of flows. This has been successful in other 
cases 0, 13 ^ • For pipe flow Barnes and Kerswell 
showed that the instabilities that arise from rotation do 
not extend to the case of pure pipe flow. Guided by 
the dominance of downstream vortices in the station- 
ary and travelling states in planar shear flows we there- 
fore adopted the following strategy, similar to the one 
used in ;14] : at low Reynolds number a volume force F 
that generates a transverse flow with translation invari- 
ant downstream vortices was added. While there is no 
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specific physical realization underlying this volume force, 
it could be achieved with differential heating and cooling 
as predicted by the symmetries of the states, or by mag- 
netic forcing in an electrically conducting fluid '245 ■ The 
amplitude of the driving was then increased beyond the 
first bifurcation, where a stable 3-d TW without transla- 
tional symmetry appears. A path following scheme was 
then used to track this forced state through parameter 
space: the Reynolds number was increased and the am- 
plitude of the body force decreased. If a situation with 
zero body force could be reached, a 3-d TW for the orig- 
inal pipe flow has been found. 

In this way we have identified four different states that 
have a discrete azimuthal rotation symmetry by construc- 
tion. An n-fold symmetry C„ is defined by the invariance 
under rotation around the pipe axis by an angle 27r/n, 
i.e. 

(ur,M0,u^)(r, 0, z) (ur,U0,Uz)(r, (/) + 27r/n, z). 

The resolution in azimuthal direction was adjusted to 
the fundamental domains of angle 27r/n to avoid carrying 
vanishing amplitudes. 

The state that bifurcates at the lowest Reynolds num- 
ber has a three-fold rotational symmetry in azimuthal di- 
rection (Fig. The dynamics of the states is governed 




FIG. 1: Travelling wave with threefold azimuthal symmetry 
at the bifurcation at Re = 1250. The frames are cross sections 
at different downstream positions separated by Az = Lz/6. 
Only half a period is shown: the last frame is the same as 
the first one up to a reflection at the horizontal diameter 
{(j> — > —(j),u^ —u^). Note that the six high speed streaks 
near the wall move much less than the low speed streaks closer 
to the center. Velocity components in the plane are indi- 
cated by arrows, the downstream component by color cod- 
ing: velocities faster than the parabolic profile are shown in 
green/yellow/red, slower ones in blue. 



by vortices which have a predominant downstream ori- 
entation and are slightly tilted inside the volume. They 
resemble near wall coherent structures in turbulent wall 
flows T^. They transport slow fluid towards the cen- 
ter and fast fluid towards the wall, thus producing high- 
speed streaks near the wall and low-speed streaks towards 
the center (all speeds are relative to the laminar profile 
for that Reynolds number). The steeper gradients near 
the wall imply higher friction losses and higher pressure 
gradients for the same flow speed. The position of the 
fast-speed streaks hardly changes in downstream direc- 
tion, whereas the low-speed streaks change shape and 
move in position quite a bit. 

In addition to the C„ symmetry externally imposed by 
the forcing the TW's dynamically develop a 'shift-and- 
reflect' symmetry, in agreement with the streak instabil- 
ity of other model flows 25] : if is the wavelength of 
the state, then 

(1) 

The same shift-and-reflect symmetry is present in 3-d co- 
herent states in plane Couette, Taylor-Couette and plane 
Poiseuille flow 13, 14, 15, iqj. 

Each TW is part of a two-parameter family of states 
that exists for a range of Reynolds numbers and a range 
of downstream wavelengths (Fig. The lowest critical 




FIG. 2: Dependence of the critical Reynolds number on the 
downstream wavenumber for the two-, three-, and fourfold 
symmetric state. The numerical values (dots) are interpolated 
by a 4th-order polynomial fit (continuous lines). 

Reynolds number is obtained for a wavelength of about 
4.2i? for the twofold symmetry and about 2.5i? for the 
three- and fourfold symmetric state. This period seems a 
bit shorter than that for optimal states in plane Couette 
flow where the optimal wavelength is about 2TTd with d 
the gap width. However, when the wedge shape of the 
boundary is taken into account and the lengths are com- 
pared to the widths, the comparison is more favorable: 
the ratio widths to lengths is about 27r:47r — 1:2 in plane 
Couette flow and about 1:2.5 in pipe flow. 
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Other states that we could identify in pipe flow in- 
clude those with twofold, fourfold and fivefold symmetry 
(Fig. 13) . All these states have a similar topology: they 





FIG. 3: Travelling waves with symmetries d with n = 2, 3, 
4, and 5. In order to highlight the topology of the states all 
states are averaged in downstream direction. The represen- 
tation of the velocity field by vectors (in-plane motion) and 
color (downstream component) is as in Fig. [H The absolute 
scale for the velocity fields is given in Table Q 

have 2n streaks of fast fluid close to the wall and streaks 
of slow fluid towards the center. The high-speed streaks 
near the wall remain fairly stationary over one period of 
the wave, and the low speed streaks in the center oscil- 
late vigorously. All states have the shift-and-reflect sym- 
metry The wavelength that gives the lowest critical 
Reynolds number for these states varies with the symme- 
try (Fig.[2Jl. Several properties of these states are listed in 
Tabled The lowest critical Reynolds number is obtained 
for the TW with C3 symmetry, for which the arrange- 
ment of vortices is optimal in the sense of being closest 
to a hexagonal packing, the preferred pattern in other 
systems. The states with C4 and C5 symmetry are not 
included in the table as they seem to be the ones which 
are most sensitive to the numerical resolution and not 
fully converged. A preliminary estimate for their critical 
Reynolds number is about 1590 for C4 and 2600 for C5. 

Interestingly, the list of TWs does not include a state 
with the Ci- symmetry of the mode that shows the 
strongest linear transient growth j23, . A state with 
Ci symmetry could be found for non-vanishing volume 
force, but up to Reynolds numbers of 4000 it could not 
be continued to pure pipe flow without reconnecting to 
the laminar profile. This should have implications for 
low dimensional model building where linear ar- 

guments have led to an emphasis of the Ci mode. 



symmetry 


C2 


C3 


Rec 


1350 


1250 




4.19 


2.58 


c/U 


1.43 


1.29 


riu 


2 


1 


Us/U 


0.38 


0.35 


ut/U 


0.035 


0.046 



TABLE I: Selected properties of travelling waves at the 
saddle-node bifurcation. Given is the critical Reynolds num- 
ber Rcc at the optimal wavelength Lc, the phase velocity c 
and the number n„ of unstable dimensions. Us is the max- 
imum deviation of the streamwise velocity from the laminar 
flow, ut is the maximum transverse velocity component. 



The phase speed c of the TWs is a function of both 
parameters, length and Reynolds number. In all cases 
the phase speed is larger than the mean speed but slower 
than the maximal speed possible with a laminar fluid 
(Fig. 0)) . That is to say, the wave still propagates down- 
stream when viewed from a frame of reference moving 
with the mean flow velocity. 
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FIG. 4: Phase velocities for the travelling states at optimal 
wavelength as a function of Reynolds number. The veloci- 
ties are normalized by the mean downstream velocity for the 
states. Three different resolutions are plotted, indicated by 
(J/K) where J is the number of independent Fourier modes, 
K is the number of Legendre polynomials. 

By advecting the unperturbed profile small transverse 
components can produce strong downstream streaks. Ta- 
ble shows that for the TWs the transverse components 
are about an order of magnitude smaller than the differ- 
ences between the downstream components of the lami- 
nar profile and the TW. 

The TWs appear as finite amplitude saddle-node bi- 
furcations. Away from the bifurcations there is an upper 
(nodal) and a lower (saddle) branch. The saddle has one 
more positive eigenvalue than the node, but in all cases 
we found that even at the point of bifurcation at least 
one eigenvalue is positive: all states are thus unstable 
right from the point of creation. This is in contrast to 
plane Couette flow, where at least the lowest state has 
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an interval of stability, although a tiny one [l 
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FIG. 5: Travelling waves and turbulent states in pipe flow. 
The friction factor as defined in Eq. @ is shown vs. Reynolds 
number. The TWs are labeled by their azimuthal symmetries. 
Three different resolutions are plotted as explained in Fig. 0| 
The lower dashed line indicates the strict lower bound from 
the laminar profile (where / = 64/ Re), the upper dashed 
line gives the mean friction factor from a logarithmic profile. 
The full circles are experimental data from 28] or from our 
numerical simulations (the value at Re = 2500 with error 
bars) . 

The friction properties of all states together with ex- 
perimental and numerical results for turbulent flows are 
shown in Fig. |S| The friction factor / is defined as 

/ = 4RAp/pU^L, , (2) 



with Ap the pressure drop over a distance along the 
pipe and p the density of the fluid. The saddle-node bi- 
furcations take place at Reynolds numbers significantly 
below the values where typical perturbations induce tur- 
bulent dynamics. At the bifurcation the friction factor is 
higher than the value extrapolated from turbulent states 
at higher Re. The upper branches in Fig.Elcorrespond to 
the lower branches in Fig. 0] states with higher friction 
have lower phase velocity and vice versa. 



Besides the TWs presented here we can expect many 
other TWs. On the assumption that narrower vortices 
will require stronger driving we anticipate that they 
will have critical Reynolds number higher than the 1250 
found for the C3 state. In addition, there are all the 
states that bifurcate from the TWs, as well as states that 
bifurcate in secondary bifurcations. Furthermore, there 
will be homoclinic and heteroclinic connections between 
these TWs as well as other periodic and aperiodic states. 
All these states together then form a chaotic saddle in 
phase space that supports the turbulent dynamics |29j |. 
One characteristic of this chaotic saddle is an exponen- 
tial distribution of life times and this is indeed found in 
numerical simulations [s^. It will be interesting to see 
to which extent these states can be used to approximate 
properties of the turbulent fiow. 

Support by the Deutsche Forschungsgemeinschaft is 
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